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ABSTRACT 

We report the initial results of a large-scale optical-near infrared survey to 
extend the known young population of the entire Orion star-forming region down 
to the substellar domain. Using deep optical I-band photometry and data from 
the 2MASS survey, we selected candidates across ~ 14.8 deg 2 in the ~ 8 Myr old 
Ori OB la subassociation and over ~ 6.7 deg 2 in the Ori OB lb subassociation 
(age ~ 3 Myr), with completeness down to O.O5M and 0.072M Q respectively. 
We obtained low resolution optical spectra for a subsample of 4 candidates in 
Ori OBla and 26 in Ori OBlb; as a result we confirmed 3 new members in 
Ori OBla, one of which is substellar, and 19 new members in Ori OBlb, out 
of which 7 are at the substellar limit and 5 are substellar. We looked into the 
presence of accretion signatures by measuring the strength of the Ha line in 
emission. Accordingly, we classified the new members as having Classical T- 
Tauri star (CTTS) or Weak Lined T Tauri star-like (WTTS) nature. We found 
that all the new members confirmed in Ori OBla are WTTSs, while 39^22% °f 
the new members in Ori OBlb exhibit CTTS-like behavior, suggestive of ongoing 
accretion from a circum(sub)stellar disk. Additionally we found that none of the 
members confirmed in OBla show near-IR color excess while 38^21% of OBlb 
members show H-K color excess. These results are consistent with recent findings 
for low mass young stars in Orion OBI. The similarity in CTTS-like properties 
and near-IR excess across the substellar boundary gives support to the idea of a 
common formation mechanism for low mass stars and at least the most massive 
brown dwarfs. Finally, we remark the discovery of two new members classified 
as CTTSs, both exhibiting W(Hcn) < —140 A, suggesting significant ongoing 
accretion. 

Subject headings: Brown Dwarf: Very Low Mass Stars: Star Forming Regions: 



1. Introduction 

One of the main goals of contemporary astrophysics is understanding the processes 
of star and planet formation. In this context, the formation of the least massive stars 
and brown dwarfs is a key issue. Because brown dwarfs (BDs) are objects with masses 



intermediate between those of s tars and planets (O.O72M > M > 0.010 M ; iBaraffe et al. 



19981 : lOppenheimer et al. Ill998l ). understanding how they form can provide important links 



between the origin of stars and planetary bodies. 
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Current observational efforts seek to establish to what extent the re is a continuity o f 
the star formation process across the substellar limit [M ~ O.O72M ; iBaraffe et al.lll998l ). 
and wether the properties of the resulting (sub) stellar populations depend on aspects like 
the surrounding environment. Observational evidence shows that for some very young popu- 
lations the kinematic and spatial distributions of very low-mass stars (VLMSs) and BDs are 
similar, that there is con tinuity of the initial mass function between VLMSs and BDs with 
masses down to 0.02M Q (ILuhman et al.ll2007l). and that BDs accret e from circum-substella r 



disks, much like low-mass stars do (e.g. IJayawardhana et al.l 120021 ; iMuzerolle et al.l 120051 ). 



These results suggest that the low mass star formation process extends across and below 
the substellar limit. However, these findings are based on observations of regions with ages 
between 1 and 2 Myr (Taurus, IC348, Trapezium and Chameleon I), much younger tha n 
the ~ 10 Myr time span in which giant planets are expected to form (ICalvet et al.l 120051 ). 
Also, these studies do not offer clear evidence on the dependency (if any) of the formation 
processes with the environment, because they compare different regions which do not share 
a common origin. 

An important and required next step is the detection, minimizing spatial biases, of 
stellar and substellar pre-main sequence populations in star forming regions that show a 
variety of environmental conditions, span an age range from ~ 1 — 10 Myr, and ideally 
share the same "genetic" pool. Such samples can offer new insights into the formation 
and evolution of BDs: a) The time dependence of indicators of disks and of disk accretion, 
over a span of 10 Myr, will offer the best view yet of how disks evolve, under a diversity 
of environments and at both sides of the substellar limit, b) The space distribution and 
kinematics of these populations will provide important constr aints to the formation models, 
e.g. the dynamical ejection scenario (IReipurth fc Clarke! 120011 ). which suggests that BDs are 
stellar embryos ejected from their birth places as a consequence of dynamical interactions 
during the early stages of their formation, resulting in them being deprived from gaining 
mass from the surrounding gas. 

The Orion star fo rming region meets al l these requirements. It is a relatively nearby OB 
association (~ 400 pc; IBriceno et al.l (120051)'). re presentative of the birthplace for the major- 
ity of stars in our Galaxy (IBriceno et al.ll2007ah. containing regions with ages ~ 1 — 10 Myr 
and differing ambient condi tions dGenzel &: Stutzkilll989l ; IBriceno et al.ll2005l ). There is also 
a molecular cloud complex (IMaddalena et al.lll9861 ). in which the you ngest star population s 
are embedded, that shares the same kinematics as the older stars (IBriceno et al.ll2007bl ). 
suggesting that the stellar populations at various ages share a common origin. Through- 
out the whole re gion there is now a signifi cant number of confirmed pre-main sequence 
stellar members (IBriceno et al.l 120051 . l2007bl ); however, substellar objects have only been 
confirmed so far in the younger, densest regions with small spatial extent, such as the a 
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2004!: ISherrv et al.ll2004!; Kenvon et al 



Ori cluster fe.g. iBeiar et al.lll999l. l200ll; iBarrado v Navascues et a 



ula Cluster (IHillenbrand fc Carpenter 



2005; 



2000 



Gonzalez-Garcia et al 



Preibisch et al 



2005; 



2003a|; ICaballero et al. 



2006D or the Orion Neb- 



Slesnick et al.l 120051 ) . In 



order to look for the faintest, lowest mass members widely distributed across the Orion 
OBI association, down to and beyond the substellar limit, we are conducting a large-scale 
optical-infrared search for BDs that spans ~ 180 deg 2 over the entire region. Our strategy is 
to combine deep, coadded, i?-band and /-band data fro m our optical photo metry obtained 
with the QUEST I CCD Mosaic Camera at Venezuela (jBriceno et al.ll2005l ). with near in- 
frared measurements from the Two Micron All Sky Survey (2MASS) to select objects, which 
we call candidates, whose photometric magnitudes and colors are consistent with those of 
BDs and VLMSs. We then obtain follow-up spectra to confirm membership and establish 
their stellar or substellar status. Finally, we characterize objects in terms of spectroscopic 
signatures such as strong Ha emission, which we use as a proxy for ongoing disk accretion. 

In this work we present initial results in this effort to extend the known young popu- 
lation of the entire Orion star- forming region throughout the substellar limit. We selected 
optical/near infrared candidates within an area of ~ 21 deg 2 , encompassing part of the Orion 
OBla and most of the Orion OBlb subassociations. We conducted spectroscopic follow-up 
of candidates in an area of 0.8 deg 2 in OBla, and 2.5 deg 2 in OBlb. The optical photometric 
observations are described in section [2j Candidate selection was done by defining regions in 
optical/near-IR color- magnitude and color-color diagrams, as explained in section [31 Section 
H] shows the results of the spectroscopic observations and spectral classification of a first sam- 
ple of candidates, covering ~3.3 deg 2 . The criteria followed to establish membership and the 
stellar or substellar nature of the candidates are presented in section [51 In section [6] we study 
the spectral types and Ha emission in order to classify the new BDs and VLMSs according 
to their Weak or Classical T-Tauri-like signatures, and compare the fraction of both types of 
objects with previo us results obtained in the stellar low-mass regime (0.25 < M/M Q < 0.9) 
of the same region (JBriceno et al.ll2005l . l2007bl ). In section 7 we discuss the infrared emission 
properties of the new members, and compare with similar mass objects confirmed as mem- 
bers in other subregions of Orion. In section [8] we estimate the extinctions, temperatures, 
luminosities and masses and plot the newly confirmed members in H-R diagrams. Finally in 
section we summarize the results and conclusions. 
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2. Photometry 
2.1. Coadding of optical images 



Since 1998 a multi-band (Johnson- Cousins BVRI) and multi-epoch large-scale survey (~ 
180 deg 2 ) of the Orion star-forming region (5h < a < 6h,—6° < 5 < 6°) is being performed 
with the Jiirgen Stock 1.0/1.5 Schmidt-type telescope and the QUEST-I camera, at the 
Venezuela Nat ional Astronomica l Observatory. The camera was developed by the QUEST 
collaboration (jBaltay et al.ll2002l ) and consists of sixteen 2048 x 2048 pixel CCDs in a 4 x 4 
array which, taking into account the gaps between CCDs, provides a field of view of 5.4 
deg 2 with a scale of ~ 1.02 arcsec/pix. The system is optimized for observations performed 
in drift-scan mode at declinations near to the celestial equator. In this mode of operation, 
the telescope is fixed while the sidereal motion is compensated by adjusting the angle of 
each row of CCDs, which ride on a common supporting Invar rod, and by fine-tuning the 
readout frequency for each chip. Each CCD row is oriented in the north-south direction and 
is provided with a filter, resulting in a quasi- simultaneous observation in four different filters 
with an integration time per filter of 140 sec at 5 = 0°. Thus, the system produces scans at 
a rate of ~ 34.5deg 2 / hour / filter . 

Our main interest is the analysis of existing /-band observations, for which the system 
attains a limiting magnitude Iu m ~ 19.5, with a completeness limit of I com ~ 18.5. In order 
to evaluate the faintest magnitude at which we could expect the substellar limit to correspond 
for each of the regions spanned by our optical data, allowing for uncertai nties in ages, distanc e 
and reddening, we made the following estimates. Using the models of iBaraffe et al.l (119981 ). 
the substellar limit for OBla on the 12.6 Myr isochrone is placed at 7~ 18 . 1 (for a distance 
modulus (m — M)oBia = 7.59 and extinction Ay = 0.6; iBriceno et al.ll2005l l2007bl ). and the 
substellar limit for Ori OBlb on the 6.3 Myr isochrone is place d at I ~ 18.1 too (distanc e 
modulus (m — M)oBib — 8.21 and maximum extinction Ay = 2; lBricefio et al.ll2005l l2007bl ). 



We consider these isochron es as rough upper l imits for the age of each subregion based on 
the age ranges reported by (jBriceno et al.ll2005l ) for OBla and OBlb. Because the substellar 



limit in each region is near the completeness limit of our individual I-band observations, we 
use a coadding technique in order to increase the signal to noise ratio (SNR) and therefore 
the limiting and completeness magnitudes of our optical data. 

For the coadd we selected 8 individual (single) scans centered at 5 = —1°, that were 
obtained in the R and / bands between December 1998 and January 1999. The corresponding 
observation log is shown in Table Olfl The sum was performed using the packages Offline 



1 This publication makes use of the Sistema de Coleccion de Datos Observacionales del Telescopio Schmidt 
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( IBaftay et al.ll2002l ) and DQ developed by the QUEST collaboration. Offline automatically 
processes an individual scan considering each CCD of the array as an independent device, 
reducing the images by bias, dark current and fl at field, using calibr ation images obtained in 
drift-scan mode during the observation runs (see lBaltay et al.l (120021 ) for further details on the 
reduction process). After these corrections, Offline performs the detection of point sources, 
aperture photometry, and finally computes the ast rometric matrix for each CCD image, 
based on the USNOA-2.0 catalog (IMonet et al.lll998l ). DQ processes the raw scans, regarding 
each CCD from different scans as independent devices; then does the bias subtraction, dark 
current and flat field corrections, following the same procedure as employed by Offline. 
Using the astrometric matrices produced by Offline, the DQ software computes the offsets, 
rotations and scale differences between images of the same area of the sky obtained from 
different scans. With these corrections DQ produces, for each scan, a new output "scan" 
composed by a set of processed images that can be added pixel by pixel to produce the final 
coadded scan. The final gain in depth was ~ 1 magnitude and the resulting coadded scan 
encompasses the region 75° < a < 85°, —2.2° < 5 < 0.1° for the R and I filters. 

Our first step in defining the photometric candidate sample is the selection in I vs I-J 
diagrams of objects located above the isochrone that defines the upper age limit for each 
subassociation, because it is in these regions of the color-magnitude diagrams that the pre- 
main sequence populations are expected to fall. Because of the differences in age and distance 
we need to apply this procedure independently for each subassociation, which requires the 
definition of a s patial limit betwee n the two. Using the integrated 13 CO emissivity observed 
in the r egion bylBallv et al. (J19871 ). and the isocontours corresponding to the dust extinction 



map of ISchlegel et al.l (119981 ). as indicators of the position of the molecular clouds associated 
with OBlb, we considered a ~ 82deg as the rough limit between OBla and OBlb within 
the 2.3 deg wide strip spanning the photometric sample analyzed in this wor k. This spatial 



l imit between the two subassociations is consistent with that discussed by iBriceno et al. 



( 120051 ). who consider the molecular ringlike structure roughly centered on the Orion Belt 
star e Ori, as a good tracer of the extent of the OBlb subassociation; within this gas ring are 
found the majority of low-mass PMS stars classified by them as members of OBlb. Figure 
[T] shows part of the entire area of our Orion optical survey, as well as the strip centered at 
<5j2ooo.o = —1° considered here, the dust extinction map and 13 CO emissivity used to define 
the limit between both regions and the fields of the spectroscopic followup analyzed in this 
work. 



version 2000 (SCDObs2000; IPonsot et alj|2007l ) 
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2.2. Optical photometry and 2MASS data 

The detection of point sources in the coadded scan was performed with the daof ind 
task in IRAFq extracting objects with R and /-band fluxes above 3a s k y - The small amount 
of crowding in our images (~ 1.4 sources/arcmin 2 ) justifies the use of aperture photometry. 
We used the phot task in IRAF, adopting an aperture radius (computed from growth curves) 
of 4", which is ~ 1.4 times the typical FWHM in both bands. Because our pixel scale is 1" 
and the typical FWHM of the images is 2.5" — 3", each stellar profile is well sampled. 



The astrometric solutions were computed using the WCStools package (JMinklll999l ). and 
astrometric standards of the Guide Star Catalog GSCIlQ yielding a median offset between 
_R-band and /-band catalogs of 0.89" with a standard deviation of 0.78". 

The photometric calibration to the Cousins system was performed by defining a set 
of 416 secondary photometric standards, located in four fields within our <5/2ooo.o = — 1° 
Orion strip. We made sure that each of the 4 strips at constant declination that compose a 
single driftscan observation, as produced by each of the 4 rows of detectors in the QUEST-I 
camera, would enclose roughly 100 of these secondary standard stars. The stars were chosen 
as objects known not to be variable from our multi-epoch survey data; only those flagged 
as non-variable at the < 0.05 mag level were used. The secondary standards span a range 
of brightness but are not saturated in any single observation, and they are not too faint, 
so that each star would yield a S/N ratio > 50. The secondary standards were in turn 
calibrated usin g 12 primary photometric standard stars located in the SA92 and PG0231 



Landolt fields (lLandoltlll992l ). observed at various airmasses under photometric conditions, 
with the 4-SHOOTER CCD camera on the 1.2 m telescope at the Smithsonian Astrophysical 
Observatory. The seeing was ~ 2" throughout the observations. The 4-SHOOTER camera 
contains four 2048 x 2048 Loral CCDs separated by 45" and arranged in a 2 x 2 grid. After 
binning 2x2 during readout, the plate scale was 0.67" pixel -1 . In order to achieve more 
uniform and consistent measurements, we placed all the standard stars on the same CCD 
detector (chip 3). Bias and flat field corrections were applied to the raw images using the 
standard tasks in the IRAF ccdproc package. We then performed aperture photometry with 
the apphot package, and derived zero points, extinction coefficients and color terms using 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

3 The Guide Star Catalog was produced at the Space Telescope Science Institute under U.S. Government 
grant. These data are based on photographic data obtained using the Oschin Schmidt Telescope on Palomar 
Mountain and the UK Schmidt Telescope. 



the photcal package. The photometric errors of the secondary photometric standard stars 
are a < 0.05 in the R and I bands. 

The calibration of sources from the coadded images was performed independently for 
each CCD using standard IRAF tasks, yielding an RMS = 0.037 mag and 0.042 mag in 
the R and I bands respectively. Sources without detections in the R band were calibrated 
independently in the / band using zero points only. Therefore, we obtained two sets of data 
calibrated in different ways: one including an R-I color term and the other without it. In the 
transformation equation used to convert instrumental magnitudes to the secondary standard 
star system, the R-I color coefficient has a rather small mean value of —0.046 for the I band. 
The observed R-I color of the candidates are in the interval 1.24 < R — I < 2.4, therefore, the 
difference between I magnitudes computed using a color term and those computed without 
it is 0.05 magnitudes for the bluest and brighter candidates (expected to have an ~ M3 
spectral type) and 0.11 magnitudes for the reddest an fainter candidates (expected to have 
an ~ Ml spectral type). The final optical data catalog includes ~ 223000 objects with 
typical instrumental errors of 0.04 < 07 < 0.2 at 19 > I > Iu m and 0.05 < ctr < 0.2 at 
20 > R > Rum. Therefore, the calibration using secondary standard stars offers adequate 
photometric accuracy for the VLMSs and substellar candidate selection and the subsequent 
analysis presented at this study. The 3<r limiting magnitudes obtained are Ru m = 21.5 
and Ium = 20.7, with completeness limits of R CO m = 20.3 and I com = 19.0. We define 
here completeness as the value at which the magnitude distribution departs from a linear 
behavior. Saturation occurs at R sat ~ I sat ~ 13.5. 

A total number of ~ 127000 sources were observed in our survey region by the 2MASS 
project, with limiting magnitudes of Ju m = 18.9, Hu m = 18, Ku m = 17.4 and completeness 
limits of J com = 16.6, H com = 16, K com = 15.4, computed following the same procedure 
used for our optical data. Typical uncertainties in the 2MASS data are 0.1 < oj^h,k < 0.5 
between completeness and limiting magnitudes for the three bands. A total number of 
~ 106000 objects from the optical catalog were identified in the 2MASS data, resulting in a 
median offset between both catalogs of 1.41" with a standard deviation of 0.89". 



3. Candidate selection 

The candidate selection was performed using optical-infrared color-magnitude and color- 
color diagrams. This method has proved successf ul in identifying VLMS and BD candidate s 



color diagrams, inis metnod nas proved successl ul in ldentilymg VLlylj and nU candidate s 
in young clusters and various star- forming regions ( Briceno et all2002uLuhman et al.ll2003al ). 



We compared the optical and infrared magnitudes and colors with isochrones and evo- 
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lutionary tracks from iBaraffe et al.l (119981 ) models, in diagrams that included /, J, H and 



i^-band data. As mentioned in section [2TTI we assum ed a distance modul us of (m— M) B\a = 
7.59 for OBla, and (m - M) OB ib = 8.21 for OBlb fariceno et al.ll2005h . The first step was 
the selection in the I vs I — J diagram of objects placed above isochrones assumed to be 
a reasonable upper age limit for each subassociation, 12.6 Myr for OBla and 6.3 Myr for 
OBlb ( §2.ip . In addition, because we are interested here in the lowest mass PMS stars and 
BDs, we requ i red th at candidates were located below the 0.3 M evolutionary track of the 
Baraffe et al.l ( 19981 ) modelgj. These diagrams are shown in Figure [2J The objects selected 
from the I vs I — J diagrams were then plotted on the H vs. I — K plane. There we checked 
their positions with respect to the same isochrones and evolutionary tracks as we used in the 
I vs I — J diagrams, in order to verify the consistency of the first selection step. Finally, in 
order to reduce the contamination from reddened early type (and therefore more massive) 
stars, we selected only objects located above the reddening line corresponding to a spectral 
type M5 in the I — K vs J — H diagrams shown in Figure [3j On average, only objects later 
than M5 (M < O.1M ), at different reddening values, are expected in this region of the 
diagram as we show in figure 3. 



The reddening lines and vectors were computed according t o the ICardelli et al.1 (119891 ) 
extinction law with Ry = 3.1 (representative of OB as sociations'). IBaraffe et al.l (119981 ) mod- 



els, and temperature-spectral type relationships from lLuhman et al.l (J2003bl ). The resulting 
catalog includes 64 VLMSs candidates in OBla and 118 in OBlb, and 58 BDs candidates 
in OBla and 77 in OBlb. Based on the completeness limits in the / and J bands, and 



Baraffe et al.l (119981 ) models, we estimate the photometric sample to be complete down to 
members with masses 0.05 M & with Ay < 0.6 for OBla, and masses 0.072 M Q for OBlb with 
Ay < 2 as we show in Figure [2j Both reddening values are representative of the typical upper 
limits measured on our spectroscopic sample (s ee below) and in low mass stars confirmed in 
both subregions by IBriceno et al.l (120051 . l2007al ) . 



4. Optical spectroscopy 

4.1. Spectroscopic observations 

Though the photometric candidate selection technique is particularly sensitive to young 
VLMSs and BDs, a fraction of ~ 25 — 30% of the sample is expected to be composed 



4 The census and study of t he Orion OBI PMS stars with masses above 0.3 M is the subject of works 
like IBriceno et all (J200lL l2005h . 
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of background and foreground objects, conta minating the ca ndidate lists. Actually, the 



field dwarfs 2MASS 05395200-0059019 (L5V; iFan et all 120001 ) and the 2MASS 05012406- 
0010452 (L4V; Reid et al. in preparation) were detected in our photometric survey with 
R = 20.27 ± 0.16, I = 17.32 ± 0.01 (R-I = 2.95) and R = 20.14 ± 0.08, I = 18.23 ± 0.02 
(R — I = 1.91) respectively and are included in our candidate list. Also, background late- 
type giants may be present. Therefore, spectroscopic confirmation of membership to the 
association is necessary. 

We obtained low resolution spectra of a first sample composed of 30 objects: 4 candidates 
located in OB la and 26 in OBlb. We use d the Hectospec Multi-fiber spectrograph on the 



6.5 MMT telescope (IFabricant et al.lll998l ). In its //5 configuration the telescope offers a 
~ 0.8 deg 2 field of view, in which 300 fibers of the spectrograph can be placed. We used 
the 270 groove mm -1 grating, that provides a spectral resolution of 6.2 A with a spectral 
coverage from 3700 to 9150 A. Each fiber subtends 1.5" on the sky and requires accurate 
coordinates for positioning. We used 2MASS coordinates for all our candidates. 

The observations were performed as part of the spectroscopic followup of the CIDA 



Variability Survey of Orion (jBriceno et al.ll2005l ) in which the main goal is the identification 



of low mass stars in the range 0.2 < M/M & < 0.9. Because of this, the spatial coverage in 
each OB la and OBlb was dictated by the number of fields observed with Hectospec in each 
region for the shallow survey, out of which only fields containing VLMSs and BDs candidates 
from our deep coadd are considered here. The 4 candidates located in OB la all fall within 
one Hectospec field, spanning ~ 0.8 deg 2 . The 26 objects in OBlb are distributed over a 
larger area of ~ 2.5 deg (three Hectospec fields). For all fields the total integration time 
was 3 x 900 sec. Because in Hectospec all objects placed in the same field have the same 
integration time, many of our faint VLMSs and BDs candidates had low S/N. However, this 
did not affect in a significant way the spectral classification as we explain in the following 
section. 

All the spectra were processed, extracted and wavelength calibrated by S. Tokarz at the 
CfA Telescope Data Center, using customized IRAF routines and scripts developed by the 
Hectospec team. A typical number of 37 fibers per field were used to obtain sky spectra; 
these were combined and the resulting spectra subtracted from the target spectra. The fields 
in the OBlb subassociation were selected far enough from the the molecular clouds to avoid 
problems due to the high sky background caused by nebulosity, however, some objects placed 
in the easternmost OBlb field do show some contamination by cloud emission lines (however, 
these lines did not affect the spectral classification). The spectra were wavelength calibrated 
using arc lamp data. The appropriate correction for the sensitivity function of the detector 
was applied, and the effects of the atmospheric extinction were computed, using standard 
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IRAF routines based on observations of the spectrophotometric standard star BD+284211. 
Light sources inside the fiber positioner of Hectospec contaminate ~ 1/3 of the fibers at the 
red edge, starting at ~ 8500 A (N. Caldwell, private communication), therefore, during the 
analysis we ignored data at wavelengths beyond this value. The observations log is shown 
in Table [3 



4.2. Spectral classification 

We performed the spectral classification based on the comparison wit h of temperature 



sensit ive sp ectral indices. We extende d the semi-automated procedures of iHernandez et al. 



( 120041 ) and ISicilia-Aguilar et all (120051 ) for their application to VLMSs and BDs, by adding 



indices measured for the 16 spectral features li sted i n Table [31 and a library of standard 



spectra of field dwarfs from iKirkpatrick et all (119991 ) . Our set of spectral indices span a 



wavelength range 4775-7940A, including spectral features such as TiO absorption bands, 
that increase in strength from K5 to M7 spectral types, and VO absorption bands, that 
start dominating spectra from M7 to early L types. 

Spectral types for noisy spectra were checked directly with the standard star spectral 
sequence using the visual extinction as a free parameter to get a more accurate result. The 
classification u is shown in Table IU and Figure @] shows the spectra, with adopted spectral 
types of the new members shown in labels. 



5. Association membership and substellar nature 

5.1. Membership diagnostics 

Candidates classified as late-M type objects can be divided into three different groups: 
first, the VLMSs and BDs that belong to the Orion OBI association; second, background 
stars composed by late type giant stars; and third, foreground, old M-dwarf field stars. 
We have used several criteria to assign Orion OBI membership to each candidate. These 
include the location above the ZAMS in color-magnitude diagrams, Ay values close to the 
characteristic absorption towards each region, and spectroscopic indicators that characterize 
young, low-mass objects. 



^This procedure has been done using the code SPTCLASS, available at 



http://www.astro.lsa.umich.edu/~hernandj/SPTclass/sptclass.html 
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First, we conside r the Ha A6563 line in emission, wit h equivalent width incr easing; toward 
later spectral types (IBarrado y Navascues et all l2003bl ; IWhite fc Basril 120031 ). Strong Ha 
emission is a common feature of chromospherically active young objects, and of young objects 



accre ting from a circumstellar disk (for the largest equivalent widths; e.g. iMuzerolle et al. 



20051 ). By adopting this criterion we eliminated contamination from background giant stars. 
Still, Ha emission can occur in dMe stars. In order to distinguish PMS objects from these 
late type field stars, and considering that PMS objects are still contracting, we chose as an 
additional youth indicator the strength of the Na I A8195 absorption feature. This doublet 
is sensitive to surface gravity and varies significantly between dwarf field stars and PMS 
objects; eq uivalent widths in you ng M dwarfs fall between those of field M dwarfs and M 
giant stars (JLuhman et al.ll2003bl ). Figure [5] shows the Na I lines of M type objects co nfirmed 
as members, superimp osed on spectra of field dwarfs from iKirkpatrick et all (119991 ). and of 
young M dwarfs from iLuhman et all (J2003bl ) of the same spectral type. All these spectra 
were smoothed to the same resolution of 16 A, and normalized in the range 8000-8400 A. 
Both criteria used simultaneously offer a good membership indication. Additionally, the Av 
values obtained in this sam ple are in good agreem ent with the characteristic values measured 
in both subassociations by lBriceno et all (j2007bl ) for low mass stars. Summarizing, M-type 
objects were classified as members if they showed Ha emission, Na I absorption clearly weaker 
than in field M dwarfs, and Ay consistent with the known values for each subassociation. 
The confirmed new members (3 in OB la and 19 in OB lb) are shown in Table HI In the 
OBla candidate sample we identified one late type (M6) object from the field, and in OBlb 
we classified 7 objects as field stars, 4 with late spectral types (M3, M5, M6 and M7) and 3 
reddened early type stars. 

Taking our entire spectroscopic sample of 30 objects, we find that 8 were classified as 
field stars (1 out of 4 in OBla and 7 out of 26 in OBlb). App lying the exact test for th e 
success rate in a binomial experiment (R-Statistical Software, llhaka fc Gentleman! Il996l ). 
we find that the proportion of objects classified correctly as members is 0.73, with a 95% 
confidence interval of 0.54-0.87. Therefore, the selection technique described here identifies 
correctly 73.3% of objects. 



5.2. Substellar status 



According to the models of iBaraffe et all (119981 ) . objects with masses at the substel- 
lar limit (0.072 M ) and ages between 7.9 Myr and 3.2 My r (that correspond to the mean 
ages of OBla and OBlb computed by iBriceno et all (120051 ) based on PMS low-mass stellar 
populations), have effective temperatures of 2936 K and 2993 K respectively. Therefore, 
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adopting the temperature to spectral type relationship from iLuhman et al.l (J2003bl ) the sub- 
stellar limit at these ages corresponds to a spectral type between M6 (2990 K) and M6.5 
(2935 K). Following this criterion, we consider as BDs those members with spectral types 
later than M6. In the Hectospec field located in the OBla subassociation we confirmed 1 
VLMS with spectral type M4.5, 1 VLMS at the substellar limit with spectral type M6, and 
1 BD with a spectral type M7. In the Ori OBlb subassociation we identified 14 VLMSs, out 
of which 7 have M6 spectral types, and 5 BDs (two with spectral types M6.5, two M7 and 
one M7.5). 



6. T Tauri signatures across the substellar limit 



The paradigm for low-mass star formation draws a picture in which a central low-mass 
star (T Tauri star) accretes mass from a circumstellar gas and dust disk, via magnetospheric 
accretion. Meanwhile the inner regions of the circumstellar disk are disrupted by the stellar 
magnetosphere, which then channels visco usly accreting material out of the disk plane and 
towards the star along magnetic field lines (JKonigl.lll99ll ; IShu et al.l 1 19941 ). Signposts of this 
process, which have been directly observed in CTTSs, include: a) infrared emission from dust 
at a range of temperatures in the disk heated by stellar irradiation and yiscous dissipation 



(from ~ 1500°K at a few stellar radii, d own to < 50° K at ~ 100 AU, e.g. Mever et al.lll997l : 



Hartmannl 11998c iMuzerolle et al.ll2003l ); b) blue/UV continuum excess emission from the 



accretion shock formed as accreting material fa lls onto the stellar surface (IValenti et al.lll993 



Hartigan et al.lll995c ICalvet fc Gullbring]|l998); c) broad permit ted emission lines produced 



in the ballistic magnetospheric gas flows (IMuzerolle et al. 



lines produced in accretion-powered winds and jets (IHartigan et al.lll995l ). 



2001) ; and d) forbidden emission 



Recent studies have searched for similar characteristics in the lowest-mass CTTSs and 
young BDs, as one way to ascertain whether similar formation mechanisms apply between 
both types of objects. There is indeed considerable evidence for magnetically-mediated disk 
accretion, including broad permitted emission line profiles (especially Ha) and optical con- 



tinuu m veiling in BDs (e.g. IWhite fc Basrill2003l ; IMuzerolle et al.l 12005c iJayawardhana et al. 



20031 ). Furthermore, the infrared exc ess emission detected in many objec t s confirms the pres- 
ence of irradiated circumst ellar disks (IMuench et aU2001clLiu et al.ll2003c I Jayawardhana et al 
20031 ; iMohanty et al.ll2004l ) and high resoluti on observations of Ha emission allowed the mea- 
surements of the mass accretion rates in BDs ( Muzerolle et al. 120001) resulting i n a much lower 



value than stellar- mass CTTSs with similar ages (e.g. IMuzerolle et al.ll2005c ILuhman et al. 
20071 ). It now seems well established that there is a continuum of CTTS-like properties 



across the substellar limit, shared by both the lowest-mass CTTSs and at least the higher 
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mass young BDs. 

Low-resolution spectra provide an appropriate means to look for Ha emission in CTTS. 
The dividing line between accreting (CTTSs) and non-accreting (WTTSs) was first set 
in terms of Ha emis sion at W(Ha) = —10 A by iHerbig fc Belli (119881 ). More recently, 
White fc Basril (120031 ) have revised this classification. In their new scheme, a star is consid- 
ered to be a CTTS if W(Ha) < -3 A in the range K0-K5, W(Ha) < -10 A for K7-M2.5, 
W(Ha) < -20 A for M3-M5.5, and W(Ha) < -40 A for M6-M7.5. Here we followed their 
criteria to classify an object as a CTTS or a WTTS. Additionally Figure 6 shows the Hel 
emission for the spectra in which this line was detected. Figure [7] shows the classification, 
with W(Ha) plotted as a function of the spectral type, with the CTTS/WTTS divide from 
White fc Basril (120031 ) indicated by the dotted line. The results are summarized in Table HI 



7. Disks among the widely spread PMS population of Orion OBI 
7.1. Trends in accretion-related indicators 



The fraction of objects exhibiting CTTS-like properties can be used as an indicator of 
the population of disk-bearing objects which are still accreting from circum(sub) stellar disks. 
How these numbers compare among samples with differing ages, but presumably sharing a 
common origin, can provide insight into the time scale for the shutting off of accretion near 
and below the substellar boundary. 

Using the same disk indicators we consider here, the first determination of the CTTS and 



WTT S fractions for the low-mass, widespread stellar population in Orion by IBriceno et al. 
(120051 ). over an area of ~ 68 deg 2 showed that 11% of the m e mbers in OBla and 23% of 



the members on OBlb are CTTS. More recently IBriceno et al.l (l2007bl ) found a lower CTTS 
fraction of 6% in Ori OBla and 13% in Ori OBlb. However, this study was restricted to a 
much smaller area (totaling ~ 1.6 deg 2 ). It is worth noting that because there are regions 
with higher concentrations of CTTS, like the very young c lusters NGC 2024/2 068 located 
on the Orion molecular clouds (e.g. see Figures 5 and 6 of IBriceno et al.ll2005l ). the CTTS 
fraction on small spatial scales is strongly dependent on the location of the region being 
considered. 

The very small size of our sample in OBla (only 3 members) does not allow us to derive 
any meaningful estimate, other than saying that it is completely dominated by WTTSs. In 
OBlb, we find that 7 of the 18 members that could be assigned a TTS type, are CTTS 
(because of suspected contamination from nebular emission, 05390910-0200267 could not be 
classified as either WTTS or CTTS). If we again apply the exact test for the success rate in 



-15- 



a binomial experiment (R-Statistical Software, llhaka &: Gentleman! Il996l ). we find that the 
proportion of CTTS in OBlb is 0.39, with a 95% confidence interval of 0.17 - 0.64; the large 
size of the confidence interval being due to the relatively small sample size. These results are 
largely consistent with previous findings, in the sense that the OB la population is dominated 
by non- accreting PMS stars, whereas the accretor fraction in the younger OBlb region can 
be estimated in the range ~ 13 — 40% depending on the region being considered and the 
sample size. 

Our findings suggest that the overall number of accretors amon g VLMSs and BDs i n 
OBlb is similar to that determined for low-mass (M > 0.3M Q ) stars by lBriceno et al.l (120051 ). 
Also, it seems that, as it happens in the low-mass PMS members, the VLMS/BD accretor 
fraction falls off by a significant amount once the population reaches ages ~ 8 — 10 Myr. 
Determining the detailed evolution of the accretion fractions is similar between low-mass 
stars and VLMSs/BDs will require larger samples. 



7.2. Near IR disk indicators 



We analyzed the J, H and K band emission of all the new members confirmed in this 
work. Figure M shows the observational J — H vs. H — K diagrams with the members 
classified as WTTSs and CTTSs. The members confirmed in OBla (age ~ 8 Myr) do not 
show excess emission in either color. In th e ~ 4 Myr old OBlb , where a fraction of the low- 
mass PMS stars (M > 0.3M Q ) identified by lBriceno et al.l (120051 ) show infrared excesses, 7/19 
members confirmed in this work (1 BD WTTS, 4 VLMSs WTTSs and 2 VLMSs CTTSs) 
have H — K > 0.45 and lie to the right of the reddening line for an M6 star in figure [HI the 
region expected for ob jects with ex c ess emission produced by hot dust in the innermost part 
of the disk (~ 0.1 AU; lMever et all Jl997l ): lMuzerolle et all (b()03h ). The resulting inner-disk 
fraction for the entire sample (VLMSs AND BDs) is 0.38, and using the exact test for th e 
success rate in binomial experiments (R-Statistical Software, llhaka fc Gentleman! ( 19961 )) 
we find a 95% confidence level interval 0.17-0.64. These numbers are very similar to what 
we found from disk accretion indicators (17.11) . and consistent with the findings reported by 
Bricefio et all (J2005J ). 



Five of the new members confirmed in this work classified as WTTS VLMSs with 
masses (0.08 < M/M Q < 0.4) are located in the y 1 deg 2 OBlb field studied with the 
Spitzer Space Telescope by Hernandez et al.1 (l2007bl ). Unfortunately, there are no similar 
observations for the objects we identified here in OBla. For the five OBlb objects we have 
photometry in the four IRAC bands with the exception of 05312373-0150245, not detected 
in the [8.0] \xm band. These five members have all been classified as WTTS and Table 5 
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provides the IRAC photometry. Figure [9] shows the [3.6]- [4.5] vs. [ 4.5]-[5.8] and [3.6]-[4.5 ] 
vs. [5.8]-[8.0] color-color diagrams with the excess region defi ned by Luhman et al. ( 20051 ) 
(left panel) and the photospheric and CTTS locus defined by lHartmann et all (120051 ) (right 
panel). All the objects exhibit photospheric colors in both diagrams, consistent with their 
optical classification as non-accreting objects; the exception is 05300324-0138428, classified 
as WTTS but whose [ 3.6]-[4.5] vs. [5.8]-[8.0] colors fall within the CTTS locus defined by 
Hartmann et al.l (120051 ). albeit offset toward smaller excess values. 



Having only four objects in OBlb (all WTTS) with full IRAC photometry prevents us 
from de riving a meaningf u l disk f raction; however, it is worth mentioning that for this same 
dataset iHernandez et al.l (J2007bl ) derived a disk frequency of 13%, using ~ 100 low-mass 
(~ Q.3M n ) PMS stars. Their determination agrees well with the CTTS fraction derived 



by iBriceno et al.l (J2007bl ) in the same area. O ur inner-disk fraction of 38_ 21 % (see above) 



lernandez et a. 



(l2007br) , is higher than their 



Briceno et al.l (120051 ). As we have discussed 



obtained in an area ~ 3 times larger that of the 

estimate, but consistent with the results from 

in section I7TTI the discrepancy could well be due, at least in part, to the largely non-uniform 

spatial distribution of disk-bearing CTTS in OBI. 



In the young a Orionis cluster (age ~ 3 Myr) ICaballero et al.l (120071 ) combined optical, 
near-IR and Spitzer IRAC/MIPS measurements to derive a disk fraction of ~ 50% in a sample 
of 30 BDs down to M ~ 0.015 M & . However, this sample is a mixture of spectroscopically 
confirmed members and candidates selected on the basis of photometry only. If we restrict 
the analysis to the 13 objects in their Table B.l which have spectral types, in the range 
M6-M7.5 (a spectral type range that overlaps the sample in this work), we find that only 5 
objects exhibit excess IR emission that can be attributed to circum(sub)stellar disks. This 
produces a somewhat lower disk fraction of 38%, which agrees very well with the value we 
found above, from the near-IR disk indicators. 

As in the case of the accretion-related indicators, our results from the JHK 2MASS 
photometry suggest that: 1) Disk fractions are similar in VLMSs/BDs and the lower mass 
PMS stars (M > O.3M ) at a given age. 2) The disk fraction near the substellar boundary 
falls off significantly, as occurs for the lower mass stars (0.8 < M/M Q < 0.3), during the 
period ~ 4 — 8 Myr. However, the present small size of the available samples precludes any 
detailed analysis on how these fractions may evolve on both sides of the substellar limit. 



We want to remark that strictly following the classification scheme by White fc Basri 



(120031 ). of the 7 objects which show IR excesses, namely 05300324-0138428, 05312373- 
0150245, 05350579-0121443, 05322069-0125511, 05340726-0149380, 05345443-0144473 and 
05390054-0150469, only the last two have been classified as CTTS. Object 05300324-0138428 
shows an Ha emission equivalent width at the WTTS/CTTS limit in Figure 7, and its IRAC 
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colors place it at the borderline between photospheric colors and IR excesses on both color- 
color diagrams of Figure 9. A similar situation is that of 05312373-0150245, which is located 
near the WTTS/CTTS boundary in Fig. 7, given our 0.5 subclass uncertainty in the spectral 
type classification, and its IRAC colors suggest it may have a small IR excess. Objects 
05322069-0125511 and 05340726-0149380 have Ha emission equivalent widths well below 
the WTTS/CTTS boundary. We have IRAC photometry only for 05322069-0125511, and 
its colors are clearly photospheric. In the near-IR JHK diagram (Fig. 8) both objects 
have the smallest H — K values among the "excess emission" objects, lying closest to the 
H - K = 0.45 limit of our JHK excess box. BD 05350579-0121443 whose J, H and K 
magnitudes are among the faintest of the sample, has an "ABB" 2MASS photometric flag 
which indicates a lower SNR in the H and K bands, and hence a larger uncertainty in the 
H — K color. We conclude that, of the 5 objects classified as WTTS but showing IR excess, 
4 have Ha equivalent widths, from our low-resolution spectra, that place them on or close 
to the border between WTTS and CTTS types. High resolution spectroscopy will be needed 
to analyze the Ha line profile and determine t he presence of line wi ngs extending to large 
velocities, characteristic of accreting stars (e.g. iMuzerolle et al.l 120051 ) . 



Of the 7 objects classified as CTTS only 2 fall within the IR excess box in Fig. 8. 
However, 3 out of the 5 CTTS outside the box have large enough 1 — a errorbars that would 
bring them into the excess emission region. Therefore, within errors, only two CTTS really 
seem to lack IR excess emission at JHK. This is consistent with fin dings for very low-mass 
stars in other star- forming regions like Taurus (JLuhman et al.ll2003al ). Observations at longer 
wavelengths should allow to detect excess IR emission in these two objects. 



8. Reddening and H-R diagrams 



The extinction towards each object was estimated from the I — J color, for which 
possible contributions from excess emission at short and long wavelengths are minimized 
(JLuhman et al.ll2003bl ). In order to compute the intrinsic colors as a function of spectral 
type, we inter polated the cor r espond ing tempe rature from the spe ctral type-temperature 
relationship of iLuhman et al.l (j2003bl ). into the iBaraffe et al.l (119981 ) models corresponding 
to ages of 7.9 Myr for OBla, and 3.2 Myr for OBlb. 

The resulting intrinsic colors and computed extinctions are listed in Table HI Extinctions 
are in the range 0.6 < Ay < 0.9 for OBla a nd 0.1 < Ay < 3.9 for OBlb consistent with 



previously reported values in both regions by IBriceno et al.l (120051 ) . Bolometric magnitudes 
were calculated using absolute / band magnitudes Mj and the bolometric correction BCi = 
0.02+0.575(l / -/)-0.155(y-/) 2 from lLegget et al. I Jl996h where V-I is the intrinsic color 
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from iKenyon fc Hartmann I (119951 ). Figure 10 shows the resulting H-R diagrams and table 4 
the mass estimates that result from the interpolation of luminosities and temperatures within 
the iBaraffe et al.1 (119981 ) models. The mean ages derived are ~ 4 ± 2Myr and ~ 8 ± 2Myr 
for O Bla and OBlb resp ectively, in good agreement with ages derived based on low-mass 
stars (jBriceno et al.ll2005l ). 



9. Summary and Conclusions 

We have presented initial results of our search for substellar objects and VLMSs in the 
widely dispersed populations of the Orion OBla and OBlb subassociations. Our results can 
be summarized as follows: 



• Using deep optical I-band photometry and data from the 2MASS survey, we con- 
structed optical-infrared color-color and color-magnitudes diagrams, in which we se- 
lected candidates across ~ 14.8 deg 2 in Ori OBI and over ~ 6.7 deg 2 in Ori OBlb with 
completeness down to O.O5M with A v < 0.6 for OBla, and O.O72M with A v < 2 
for OBlb. We obtained low resolution optical spectra for a subsample of 4 candidates 
over ~ 0.8 deg 2 of OBla, and 26 candidates within ~ 2.5 deg 2 of OBlb. 

• Through spectral signatures we confirmed 3 new members in OBla, one of which 
is clearly substellar (spectral type M7, with an inferred mass ~ 0.06M Q ) and 2 
are VLMSs, one of which is at the substellar limit (spectral type M6, with mass 



M 



0.072M m ). In OBlb we found 19 new members: 14 are VLMSs, of which 7 



are at the substellar limit, and 5 are substellar members with spectral types between 
M6.5 and M7.5. From the spectroscopy we find that our photometric candidate se- 
lection technique is highly efficient (73%) in picking out young, very low-mass and 
subst ellar members of th e Orion OBI association. Masses were estimated according to 
the (IBaraffe et al.lll998l ) models, and the less massive member (M7.5 spectral type), 
has an estimated mass M ~ O.O4M . 

T he new member s have been classified as CTTSs or WTTSs following the scheme 
of IWhite fc Basril (120031 ) . We found that all three members confirmed in OBla are 



WTTSs, while 39^ 22 % °f the members in OBlb are CTTSs. These results are largely 
consistent with recent findings, in the sense that the OBla population is dominated 
by non-accreting PMS stars, whereas the accretor fraction in the younger OBlb re- 
gion can be estimated in the range ~ 13 — 40% depending on the exact location and 
area being considered, and the sample size. Our findings indicate that in OBlb, the 
number of accretors in VLMSs/BDs is similar to that derived for low-mass PMS stars 
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(JBriceno et all (120051 )); also, that the overall number of accretors, both in low-mass 



stars and among VLMSs/BDs falls off by a significant amount by ages ~ 8 — 10 Myr. 

Of the 19 newly confirmed members of Ori OBlb, 7 (1 BD WTTS, 4 VLMSs WTTSs 
and 2 VLMSs CTTSs) exhibit excess emission in the 2MASS H-K color as would be 
expected of thermal emission from hot dust in the innermost part of a circumstellar 
disk. None of the 3 Ori OBla members show near-IR excesses. We derive an inner- 
disk fraction of 38l2i%; which is in excellent agreement with our result for the accretor 
fraction. Our inner-disk estimate for VLMSs and BDs is also consistent with the 



results from iBriceno et all (120051 ) for higher mass stellar obje cts (M > 0.3M&), an d 



with findings for BDs spectroscopically confirmed in a Ori by lCaballero et all (120071 ). 
The few Ori OBlb members with Spitzer data are all WTTS, by optical indicators; 
the IRAC photometry is largely consistent with this classification. Only one object 
shows marginal indication of excess emission at 8/im. 

As in the case of the accretion-related indicators, the near-IR indicators suggest that: 
1) At a given age, the inner-disk fraction of objects at/near the substellar boundary is 
similar to that determined in the lower-mass (M > O.3M ) PMS stars. 2) The inner- 
disk fraction seems to fall off significantly, both in VLMSs/BDs and low-mass stars, 
once the population has aged to ~ 8 Myr, as in Ori OBla. However, the present small 
size of the available samples precludes any detailed analysis on how these fractions may 
evolve on both sides of the substellar limit. 

Finally we remark the discovery of two new members 05335219-0156398 (M7) and 
05390532-0135327 (M6) classified as CTTSs, without contamination by emission lines 
from the cloud, both exhibiting strong Ha emission, with W(Ha) < —140 A, suggesting 
significant ongoing accretion. 
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RA J2000.0 (DEG) 

Fig. 1. — Spatial location of the newly discovered VLMSs and BDs. Triangles and stars 
indicate BDs and VLMSs respectively, spectroscopically co nfirmed as members within the 
Hectospec fields (solid line circles). Members classified by iBriceno et al.l (120051 . l2007bl ) as 
CTTSs are indicated with small empty squares and WTTSs with filled dots. The horizontal 
dot-dashed lines indicate the north and south limits of our scans centered at £j2ooo.o — —1-0°. 
The dashed line region marks the boundary of the OBlb and OBla regions as defined by 
Warren fc Hesserl (119771 ). The dotted vertical line indicates the adopted spatial limit between 
OB fa and OBlb subregions into the strip of our photometric sample as we e xplained in 
sectio n 12.11 The gray scale map shows the integrated 13 CO em issivitv from |Ballv et al. 
( 19871 ). The isocontours correspond to the dust extinction map of ISchlegel et al.l (119981 ). in 
E{B — V) steps of 0.1 to 1.5 mag. 
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Table 1. Observation log of the individual scans used to produce the coadded scan 



Date 


Scan 


tti(J2000) 


a/(J2000) 


Filters 


Mean Seeing (") 


1998-12-13 


504 


04:10:00 


05:51:00 


RBIV 


2.81 


1998-12-13 


505 


04:10:00 


05:51:00 


RBIV 


2.95 


1998-12-29 


501 


04:10:00 


05:51:00 


RBIV 


2.91 


1999-01-09 


527 


03:55:00 


05:53:00 


RBIV 


2.93 


1999-01-09 


528 


03:55:00 


05:53:00 


RBIV 


2.92 


1999-01-10 


529 


03:55:00 


05:53:00 


RBIV 


2.85 


1999-01-09 


528 


03:55:00 


05:53:00 


RBIV 


2.92 


1999-01-10 


529 


03:55:00 


05:53:00 


RBIV 


2.85 


1999-01-10 


530 


03:55:00 


05:53:00 


RBIV 


2.95 


1999-01-22 


501 


03:55:00 


05:53:00 


RBIV 


3.06 
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Fig. 2. — First selection of candidates. Dots represent field sources loca ted under the 
i sochr one that defines the upper age limit of each region, based on results from lBriceno et al. 
J2005J ) (12.6 Myr for OBla and 6.3 Myr for OBlb; dashed lines), or PMS candidates located 
above the O.3M evolutionary track that were rejected as candidates. Sources above the 
isochrones, between the 0.3M Q and 0.072M Q evolutionary tracks (open triangles), and be- 
low the substellar limit (0.072M Q evolutionary track; open circles), were selected as initial 
candidates. The thin solid lines represent the completeness of the photometric sample, and 
dotted lines show the completeness assuming Ay = 0.6 for OBla, and Ay = 2 for OBlb. 
These extinctions correspond to the representative upper values measured in our spectro- 
scopic sample for each region, resulting in a m ass completene s s of ~ O.O5M for OBla and 



0.072M m for OBlb. We used models from 



'o 



Baraffe et al. 



(119981 ) and a distance modu- 



lus of 7.59 for OBla , and 8.21 for Oblb (JBricefio et al.ll2005l ). The reddening vectors were 
computed using the ICardelli et al.l (119891) extinction law with Ry = 3.1, and temperature- 
spectral type relationship from lLuhman et al.l (J2003al ). 
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Fig. 3. — Final step in the selection of candidates. The symbols are the same as in Figure [2J 
Solid lines represent the reddening line for M5 objects and dashed lines the redd ening lines 
for Ml (lower) and M9 (upper) objects , computed using the lCardelli et al.l (Il989l ) extinction 
law with R v = 3.1.lBaraffe et al.l (119981 ) models, and temperature-spectral type relationships 
from iLuhman et al.l (]2003al ). On average, sources falling above the M5 reddening line are 
later than M5 at different Ay values and were selected as final candidates. 
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Table 2. Hectospec observation log. 



Date Subregion Field ID a(J2000) 5(J2000) Integration [seg] 



Apr/09/2005 


OBla 


decml60_04 


05:17:59.67 


-01:35:19.4 


2700 


Nov/03/2004 


OBlb 


decml60_07 


05:30:00.47 


-01:35:33.7 


2700 


Mar/04/2005 


OBlb 


decml60_08 


05:33:55.64 


-01:34:29.8 


2700 


Nov/19/2004 


OBlb 


decml60_09 


05:37:55.32 


-01:34:33.0 


2700 
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Table 3. Spectral features from iKirkpatrick et al.l (119991 ) used for the semi-automated 

spectral classification. 



Feature ID A c (A) AA (A) 



TiO-1 


4775 


30 


TiO-2 


4975 


50 


TiO-3 


5225 


75 


TiO-4 


5475 


75 


TiO-5 


5600 


50 


TiO-6 


5950 


75 


TiO-7 


6255 


50 


TiO-8 


6800 


50 


TiO-9 


7100 


50 


TiO-10 


7150 


50 


VO-1 


7460 


100 


VO-2 


7840 


100 


VO-3 


7940 


100 
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Fig. 4. — New members confirmed in this work. Each spectrum shows the last four characters 
of its 2MASS ID, and the spectral type followed by its classification as CTTSs or WTTSs. 
All spectra were normalized at 7500 A and are shown at their original resolution of 6 A. 
Member 0267 does not have a clear classification as CTTS or WTTS, as is explained in 
subsection 6.1. 
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Fig. 5. — The Na I A8 195 absorption lines i n the n ewly confirmed members (solid line), for 



field dwarfs (dotted line; iKirkpatrick et al. 



(1999J)), and f or young VLMSs and BDs from 



Chamaeleon I region (~ 2Myr; dashed line; iLuhmanl (120041 ) ). New members and comparison 
templates have the same spectral types. All spectra were smoothed to a resolution of 16 A 
and normalized in the interval 8000-8400A. 



Table 4. Catalog of very low-mass stars and brown dwarfs confirmed as new members of Orion OBI. 



2MASS-ID R c ± tR C Ic ± cl c 

0115233 19.225±0.047 17.121±0.010 

0153232 20.204±0.100 17.767±0.015 

0144501 18.706±0.023 16.882±0.007 

0138428 16.272±0.013 

0136587 19.049±0.037 16.989±0.008 

0150245 19.967±0.066 17.839±0.016 

0125511 18.376±0.015 16.569±0.005 

0127204 18.059±0.018 

0136234 20.081±0.071 17.819±0.014 

0112519 18.158±0.024 16.240±0.009 

0156398 18.102±0.024 

0149380 20.148±0.100 17.876±0.017 

0144473 17.020±0.005 15.773±0.004 

0121443 18.928±0.040 

0129204 20.344±0.090 18.173±0.017 

0129140 18.676±0.020 16.484±0.005 

0124270 18.617±0.026 16.701±0.007 

0155519 18.295±0.031 

0147405 18.406±0.020 16.249±0.004 

0150469 19.613±0.065 17.732±0.018 

0135327 19.604±0.041 17.619±0.011 

0200267 19.102±0.086 



W[Ha] Linos WTTS/CTTS Av Log(L/L Q ) Teff M(M ) 



05164635- 
05185578- 
05193522- 
05300324- 
05311727- 
05312373- 
05322069- 
05323839- 
05332365- 
05333251- 
05335219- 
05340726- 
05345443- 
05350579- 
05370790- 
05372198- 
05373853- 
05374145- 
05380232- 
05390054- 
05390532- 
05390910- 



4.946 


14.415 


14.127 


M6.0±0.5 


MB 


-26.8 


5.364 


14.640 


14.260 


M7.0±1.0 


M7 


-18.5 


4.864 


14.330 


14.020 


M4.7±0.5 


M4.5 


-13.3 


3.419 


12.515 


11.895 


M2.7±0.5 


M3 


-20.0 


4.944 


14.432 


14.078 


M6.0±0.5 


M6 


-16.8 


5.643 


15.104 


14.485 


M6.2±0.5 


M6 


-24.4 


4.656 


14.205 


13.709 


M4.8±0.5 


M5 


-8.8 


5.800 


15.298 


14.986 


M 5. 5 ±0.5 


M5.5 


-26.0 


5.379 


14.856 


14.458 


M6.0±1.0 


M6 


-22.3 


4.013 


13.347 


12.921 


M5.()±1.0 


M5 


-31.1 


5.736 


14.967 


14.579 


M7.0±1.5 


117 


-207.4 


5.641 


15.271 


14.766 


M5.5±1.0 


M5.5 


-10.4 


3.381 


12.511 


11.901 


M3.0±0.5 


M3 


-32.8 


6.341 


15.835 


15.229 


M6.5±1.0 


M6.5 


-18.1 


5.869 


15.564 


15.216 


M4.8±0.5 


M5 


-4.7 


4.209 


13.703 


13.333 


M6.1±0.5 


M6 


-26.8 


4.592 


14.037 


13.625 


M6.2±0.5 


M6 


-32.4 


5.987 


15.484 


15.115 


M7.0±0.5 


M7 


-47.0 


4.088 


13.487 


13.103 


M6.5±0.8 


M6.5 


-16.5 


5.476 


15.020 


14.539 


M6.2±1.0 


M6 


-67.1 


5.356 


14.784 


14.454 


M5.8±0.5 


M6 


-140.3 


6.166 


15.559 


15.413 


M7.7±l.() 


M7.5 


-93.5 



w 

w 
w 
w 
w 
w 
w 
c 
w 
c 
c 
w 
c 
w 
w 
w 
w 
c 
w 
c 
c 



0.6 


-1.81 


0.5 


-1.80 


0.9 


-1.79 


3.9 


-0.64 


0.1 


-1.63 


0.6 


-1.85 


0.2 


-1.56 


1.0 


-1.91 


1.3 


-1.67 


1.2 


-1.19 


0.3 


-1.73 


1.0 


-1.83 


2.5 


-0.78 


1.4 


-2.00 


1.4 


-1.91 


0.8 


-1.26 


0.3 


-1.47 


0.2 


-1.83 


0.1 


-1.24 


0.7 


-1.78 


0.8 


-1.72 


1.6 


-1.90 



2990 


0.07 


2880 


0.06 


3197 


0.15 


3415 


0.40 


2990 


0.08 


2990 


0.07 


3125 


0.13 


3057 


0.09 


2990 


0.08 


3125 


0.17 


2880 


0.06 


3057 


0.09 


3415 


0.40 


2935 


0.07 


3125 


0.13 


2990 


0.10 


2990 


0.09 


2880 


0.06 


2935 


0.08 


2990 


0.07 


2990 


0.07 


2795 


0.04 
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Fig. 6.— Hel emission lines observed in 05345443-0144473, 05390532-0135327 and 05374145- 
0155519 supporting their classification as CTTS. 
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Fig. 7. — W(Ha) vs. spectral type for the new members. Solid circles and triangles indi- 
cate BDs and VLMSs respectively, confirmed as new members. Squares represent members 
classified as WTTSs and pentagons indicate members class ified as CTTSs. The d otted line 
indicates the separation between CTTSs and WTTSs from lWhite &: Basril (120031 ). The cir- 
cled dot indicates member 05390910-0200267, that could not be classified as either WTTS 
or CTTS as explained in subsection 7.1. 
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Fig. 8. — Observational J — H vs. H — K diagram of the OBla and OBlb members 
confirmed at this work. Symb ols are as in F i gure [71 The dotted lines represent the 7.9 Myr 
and 3.2 Myr i sochr ones from iBaraffe et al.l (119981 ). that correspond to the ages found by 
Briceno et al.l (120051 ) on the basis of low mass stars in each sub association. The dashed line 
represents the reddening line for an M6 spectral type, and the dotted box enc loses objects 
showin g; H-K excess. The solid line is an extrapolation of the CTTS locus from iMeyer et al. 
(119971 ) for early M stars to an M6 spectral type. 
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Table 5. IRAC photometry of the new VLM members of OB lb 



ID-2MASS 


STad 


WHa 


WTTS/CTTS 




[3.6] 


ct[3.6] 


[4.5] 


[a4.5] 


[5.8] 


[<r5.8] 


[8.0] 


[o-8.0] 


05300324-0138428 


M3 


-20.0 


W 






11.180 


0.003 


10.922 


0.003 


10.785 


0.008 


10.273 


0.008 


05311727-0136587 


M6 


-16.8 


w 






13.708 


0.009 


13.635 


0.012 


13.684 


0.054 


14.297 


0.239 


05312373-0150245 


M6 


-24.4 


w 






14.464 


0.012 


14.223 


0.017 


13.995 


0.059 






05322069-0125511 


M5 


-8.8 


w 






13.528 


0.008 


13.458 


0.011 


13.501 


0.048 


13.523 


0.090 


05323839-0127204 


M5.5 


-26.0 


w 






14.611 


0.014 


14.570 


0.022 


14.533 


0.100 


15.197 


0.420 


Note: Photometry 


is from Hernandez 


et al.l J2007bll. 


Sp 


cctral typ 


es, W(Hq 


) and WTTS/CTTS 


i classification arc 


from our 


Table [4] 
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Fig. 9. — Color-color diagrams of the five OBlb members with available IRAC photome- 
try. The dashed lines in the left panel indicate the region of excess IR emission defined by 
Luhman et al. I (J2005I ) . Dashed boxes in the rig ht panel indicate the po sition of the photo- 
spheric locus and the Class II source locus from iHartmann et al.l (120051 ) . Symbols are as in 
Figure El 
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Fig. 10. — H-R diagrams of confirmed members located in Hectospec OBl a and OBlb 



fields. Symbols are as in Figure [3 Evolutionary tracks and isochrones are from lBaraffe et al. 



(J1998I ) where the 3.2 Myr and 7.9 Myr isochrones correspond to the ages of OBla and of 
OBlb respect ively, computed from the low mass members (M > 0.3M Q ) in color-magnitude 
diagrams by (JBriceno et al.ll2005l ). 
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